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ABSTRACT 


A  mathematical  model  is  developed  for  dispersion  of  a  heated  effluent 
from  an  ocean  outfall.  Input  parameters  include  atmospheric  and  oceanic 
conditions  and  discharge  characteristics.  The  model  solves  the  steady- 
state,  two-dimensional  differential  equation  for  non-conservative 
diffusion  of  heat  in  a  moving  fluid.  The  solution  is  calibrated  and 
verified  using  data  from  surveys  conducted  at  the  Southern  California 
Edison  Company  power  plant  at  Huntington  Beach,  California.  Temperature 
fields  predicted  by  the  model  are  compared  with  the  actual  fields  for 
seven  different  surveys.  These  comparisons  indicate  that  the  model  can 
be  used  to  predict  the  large  scale  influence  of  the  outfall  on  the  local 
ocean  environment. 
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I.  INTRODUCTION '  . 


A.  GENERAL 

The  existence  of  a  low  temperature  heat  sink  for  cooling  water  is  of 
extreme  importance  to  the  efficiency  of  steam  power  plants.  Since  the 
ocean  is,  in  effect,  an  infinite  fluid  heat  sink,  power  companies  find 
it  desirable  to  locate  their  plants  on  the  coast  to  utilize  the  ocean 
waters  as  condenser  coolant;. 

Steam  electric  power  plants  circulate  large  quantities  of  cooling 


water  through  their  condensers  to  extract  excess  waste  heat  from  the 


steam  leaving  the  low  pressure  side  of  their  turbines.  The  quantity  of 


waste  heat  produced  is  a  function  of  the  size  and  efficiency  of  the 
power  plants.  The  size  and  number  of  these  plants  are  increasing, 
particularly  those  using  nUclear  fuel.  Nuclear  pcwer  plants  are  even 
less  efficient  in  their  utilization  of  heat  than  fossil  fueled  power 
plants.  This  is  because  technological  difficulties  make  it  impractical, 
uneconomical,  or  unsafe  to  produce  high-pressure,  superheated  steam  in 
a  water-cooled  reactor  system  which  would  allow  more  efficient  opera- 
j  tion  [Ref.  25].  . 

Due  to  the  large  heat  capacity  of  the  ocean,  the  mean  ocean  tempera¬ 
ture  will  not  be  significantly  affected  by  the  introduction  of  this 

j 

waste  heat.  Local  temperatures,  however,  may  be  greatly  affected  and 

care  must  be  taken  to  minimize  any  adverse  effects  on  the  environment. 

j 

Hence  it  is  highly  desirable  to  develop  a  model  that  adequately  predicts 

■ 

the  temperature  distribution  and  the  area  of  influence  of  the  waste  heat. 
The  results  of  such  a  model  could  enable  estimates  to  be  made  of  the 


i 


influence  of  the  outfall  on  marine  flora  and  fauna,  and  also  aid  in 
planning  the  location  of  the  cooling  water  intake  so  that  recirculation 
can  be  kept  to  a  minimum. 

B.  EFFECTS  OF  THERMAL  POLLUTION 

Ecologists  consider  temperature  the  primary  control  of  life  on 
earth.  Each  aquatic  species  becomes  adapted  to  the  seasonal  variations 
in  the  temperature  of  the  water  in  which  it  lives  but  it  cannot  adjust 
to  the  shock  of  an  abnormally  abrupt  change  [Ref.  8].  In  addition  to 
direct  mortality  from  a  temperature  change,  there  may  be  indirect 
changes.  For  example,  temperatures  which  are  not  lethal  themselves  may 
affect  metabolism,  reproduction,  longevity  and  growth,  as  well  as  in¬ 
creasing  or  decreasing  important  food  organisms.  For  these  reasons 
there  is  growing  concern  among  ecologists  about  the  heating  of  aquatic 
habitats  by  man's  activities. 

In  an  ocean  outfall  several  detrimental  effects  could  occur.  The 
ecosystem  could  be  upset  by  hot  waste  water  to  the  extent  that  a  vital 
link  in  the  food  chain  could  disappear,  increase  in  abundance,  or  be 
replaced  by  another  organism.  If  this  occurred,  there  could  be  an 
increase  or  decrease  in  game  fish  and/or  an  increase  in  other  less 
desirable  species.  Historically,  it  is  known  that  in  the  summer, 
Macrocystis  (giant  kelp)  plants  in  the  San  Diego  area  suffer  from  pro¬ 
longed  exposure  to  warm  surface  waters  and  as  a  result  "die-back".  They 
regenerate  quickly  with  the  arrival  of  cooler  surface  water  in  the  fall 
[Ref.  17].  Not  only  is  kelp  harvested  commercially,  but  it  is  a  haven 
for  many  marine  species.  The  influence  of  a  hot  effluent  could  prevent 
the  regeneration  of  kelp  and  its  disappearance  could  completely  upset  the 
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ecology  in  the  area  of  the  outfall.  There  is  an  additional  problem 
where  a  cooling  water  outfall  is  in  the  vicinity  of  a  sewage  outfall. 
Oxygen  and  sunlight  are  important  in  neutralizing  sewage.  The  presence 
of  warmer  water  not  only  limits  the  amount  of  dissolved  oxygen,  but  also 
the  warm  water  may  be  less  dense  and  override  the  sewage  cutting  off  its 
access  to  the  atmosphere.  . . 

Beneficial  uses  of  a  heated  effluent  are  also  possible  and  should  be 
more  fully  explored.  Presently,  there  are  desalinization  plants  built 
or  planned  in  conjunction  with  power  plants  that  use  the  waste  heat  in 
their  operation.  An  increase  in  temperature  along  a  section  of  beach 
may  make  it  more  desirable  as  a  recreation  area.  Longshore  currents 
inside  the  surf  zone  could  carry  the  heat  excess  over  a  long  section  of 
coastline.  Coolant  water  could  further  be  used  on  a  varying  basis  to 
control  the  temperature  in  a  limited  area  to  allow  fish  farming  or 
aquaculture . 

C.  BACKGROUND 

A  model  to  predict  the  dispersion  of  a  heated  effluent  from  an 
ocean  outfall  should  include  all  the  physical  and  thermal  processes 
that  actually  act  on  the  effluent.  Among  the  processes  and!  (features 
that  should  be  included  are  the  volume  and  temperature  of  tins  dis- 
charge;  the  size,  depth  and  orientation  of  the  discharge  opening;  the 
thermal  and  saline  structure  of  the  ambient  receiving  water;  the 
currents;  the  waves;  the  exchange  of  energy  across  the  air-sea  inter¬ 
face;  horizontal  and  vertical  diffusion;  and  possibly  topography. 

The  heated  coolant  water,  being  less  dense  than  the  ambient  sea 
water,  rises  with  considerable  momentum.  The  column  of  heated  water 
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entrains  cooler  water  and  is  diluted  as  it  rises  [Ref.  1],  If  the 
effluent  reaches  the  surface,  it  spreads  out  as  a  thin,  warm  layer 
[Ref.  233.  The  strong  temperature  gradient  between  the  thin  layer  of 
warm  water  and  the  cooler  receiving  water  inhibits  vertical  mixing,  and 
dispersion  of  the  warm  effluent  is  primarily  horizontal.  In  the  presence 
of  a  current  the  surface  isotherms  are  roughly  elliptical;  the  ratio  of 
major  to  minor  axis  being  a  function  of  current  speed  -  the  patch  being 
elongated  in  the  direction  of  flow.  Usually  there  is  quite  a  bit  of 
irregularity  in  the  observed  isotherms,  which  may  be  attributed  to  a 
number  of  phenoma  including  random  current  motion,  internal  waves,  and 
upwelling  or  downwelling. 

In  a  classical  study  of  outfall  polluting,  Rawn  and  Palmer  [Ref.  18] 
discharged  jets  of  colored  water  into  sea  water  in  the  field  through 
nozzles  of  various  sizes  and  orientations.  From  these  experiments,  they 
were  able  to  establish  empirical  equations  relating  field  thickness  to 
the  depth  at  which  the  outfall  is  located,  the  radius  of  the  top  of  the 
jet  to  both  the  depth  of  the  outfall  and  the  diameter  of  the  outlet,  and 
dilution  to  the  distance  from  the  outfall.  Edinger  and  Geyer  [Ref.  11], 
employing  temperature  data  collected  at  several  locations  in  an  embay- 
ment  over  a  five  month  period,  obtained  a  solution  for  temperature 
distribution  which  involved  the  following  parameters:  plant  pumping 
rate,  exchange  of  heat  with  the  atmosphere,  and  apparent  diffusivity 
in  both  upstream  and  downstream  directions.  Their  solution  neglects 
any  ambient  current  and  is  based  on  the  assumption  that  the  embayment 
waters  are  uniformly  mixed  in  the  vertical  and  cross  channel  directions. 
Wada  [Refs.  27  and  28]  describes  the  effect  of  various  factors  on  the 
recirculation  of  cooling  water  in  a  bay.  He  also  presents  a  theory  from 
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which  quantitative  prediction  of  water  temperature  and  velocity  fields  may 
be  obtained  for  a  bay  using  numerical  techniques.  Pritchard  and  Carter 
[Ref.  6]  predicted  the  temperature  distribution  from  a  proposed  river 
outfall  by  using  a  fluorescent  dye  tracer  to  analyze  the  physical  pro¬ 
cesses  of  advection  and  diffusion.  A  correction  for  the  non-conservative 
nature  of  heat  was  applied  in  the  solution.  The  buoyancy  of  the  heated 
effluent  was  treated  as  an  unknown  and  solutions  were  presented  for 
several  different  depths  of  influence  of  heat. 

Silberman  and  Stefan  [Ref.  21]  used  a  hydraulic  model  to  predict  the 
effects  of  an  outfall  on  water  temperature  in  a  lake.  In  another  lake 
study  Hoppes,  Zeller  and  Rohlich  [Ref.  26]  conducted  field  studies 
throughout  an  entire  year  to  assist  in  their  development  of  a  model  to 
predict  heat  dissipation  and  induced  circulation  from  surface  outfalls. 
Other  noteworthy  research  in  these  areas  is  contained  in  References  [5], 
[9],  [10],  [13],  and  [15]. 

In  this  study,  a  mathematical  model  has  been  developed  for  heat 
dispersion  from  an  outfall  in  an  open  coastal  area.  The  techniques  and 
conclusions  are  designed  to  have  general  application  but  are  specifi¬ 
cally  based  on  part  of  a  continuing  survey  that  was  conducted  by  Marine 
Advisors,  an  oceanographic  consulting  firm,  for  Southern  California 
Edison  Company,  in  the  waters  adjacent  to  the  power  plant  located  at 
Huntington  Beach,  California  [Refs.  16  and  17].  This  survey  includes 
information  on  physical  oceanographic  conditions,  marine  flora  and 
fauna,  fish  catches,  bottom  sediments  and  suspended  fine  sediments  as 
well  as  temperature  distribution.' 

The  Huntington  Beach  power  plant  withdraws  and  discharges  its  con¬ 
denser  cooling  water  through  two  large  pipes,  each  terminating  in  the 
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/ 

ocean  at  depths  of  about  34  feet  and  28  feet  respectively.  They  are 
located  approximately  1500  feet  from  the  beach.  The  slope  of  the  bottom 
in  the  area  of  the  outfall  is  about  1/50. 
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II.  HEAT  BUDGET 


A.  GENERAL 

The  distribution  of  temperature  in  a  body  of  water  resulting  from  a 
discharge  of  warm  effluent  is  influenced  by  the  characteristics  of  the 
effluent  itself,  and  by  the  external  mechanisms  which  control  the 
exchange  of  heat  with  the  surroundings  or  environment.  The  air-sea 
interface  processes  are  jointly  referred  to  as  the  heat  budget  which 
can  be  expressed  as : 

Qs  -  Qr  -  Qb  -  Qh  -  Qe -  Q  (2-1) 

where  Qs  is  the  incident  flux  of  short-wave  radiant  energy  from  the  sun 
through  the  water  surface,  Qr  is  solar  radiation  reflected  at  the  sea 
surface,  Qb  is  the  net  energy  lost  by  the  body  of  water  through  the 
exchange  of  long-wave  radiation  between  the  atmosphere  and  the  body  of 
water,  Qe  is  the  energy  loss  associated  with  evaporation,  Qh  is  the 
exchange  between  the  body  of  water  and  the  atmosphere  in  the  form  of 
sensible  heat,  and  Q  is  the  net  gain  or  loss  of  heat  in  the  body  of 
water.  Conduction  of  energy  through  the  bottom,  heating  due  to  chemi¬ 
cal  and  biological  processes,  and  the  transformation  of  kinetic  energy 
into  thermal  energy  are  generally  neglected  because  of  their  small 
magnitude . 

In  the  present  study  an  analytical  solution  is  sought  and  hence  it 
was  necessary  to  express  the  total  heat  energy  exchange  across  the  air- 
sea  interface  as  a  linear  function  of  sea  surface  temperature  for  a 
given  set  of  atmospheric  conditions: 

Q  -  A  -  BTw.  (2-2) 
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In  order  to  adhere  to  this  restriction  the  individual  components  of  the 


heat  budget  were  linearized.  Each  mechanism  and  the  procedure  followed 
to  compute  its  contribution  are  considered  separktely.  The  results  are 
summed  to  get  the  total  exchange. 


B.  SOLAR  RADIATION  j 

The  solar  radiation  (Qs)  is  a  function  of  the  sun's  altitude  and 


the  cloud  cover.  Laevastu  [Ref.  14]  presents  several  empirical  formulae 
for  estimating  insolation  including  nomographic  aids  for  their  evalu- 

i  | 

ation.  Reflected  radiation  (Qr)  is  primarily  a  function  of  the  sun's 


altitude  and  cloudiness  although  the  state  of  the  sea  probably  plays  a 


minor  part.  Laevajsj:ti  treats  Qr  as  a  function  of  Qs .  James  [Ref.  3] 
combines  Qs  and  Qr  in  a  term  he  calls  effective  insolation  (Qs  -  Qr) 
and  presents  nomograms  for  its  determination.  Effective  insolation 
was  computed  using  Figure  1  and  Tables  II  and  III  from  James  [Ref.  3]. 
The  latitude  at  the  data  site  is  34°  North.  The  values  for  the  sun's 
declination  and  duration  of  sunlight  were  interpolated  from  Table  II. 

A  value  $  is  defined  as  latitude  minus  solar  declination.  The  cloud 
category  is  based  on  Lumb's  Cloud  Classification  (Table  III)  and  the 
cloud  conditions  given  in  Refs.  [7]  and  [24].  The  effective  insolation 

for  the  day  was  obtained  by  entering  Figure  1  with  the  value  of  $,  pro- 

|  j 

ceeding  to  the  appropriate  cloud  category,  then  vertically  to  the 

proper  number  of  hours  of  daylight  and  finally  moving  horizontally  to 

(  -  ... 

the  right.  The  effective  insolation  per  hour,  for  purposes  of  this 
study,  was  taken  to  be  1/24  of  this  daily  value.  Although  this  average 

value  is  not  a  true  value  for  most  of  the  data  periods,  it  was  felt  that 

| 

because  of  lag  time  in  the  effects  of  insolation  it  was  an  adequate  value 
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Since  this  value  is  a  constant  with  respect  to  sea  surface  temperature  it 
easily  fits  the  linear  restriction. 

C.  EFFECTIVE  BACK  RADIATION 

The  effective  back  radiation  (Qb)  can  be  estimated  by  treating  the 
water  surface  as  a  gray-body  emitter  where  it  is  some  function  of  the 
fourth  power  of  the  sea  surface  temperature.  One  of  the  more  popular 
formulae  used  is  that  developed  by  Anderson  [Ref.  2}. 

Qb  -  <4.75  x  10'9)  Tw4  (l-a+befl)  (2-3) 

where  a  -  0.74  +  0.025Ce'0,0584h 
b  -  0. 0049-0. 00054Ce'°*060h 
efl  ■  Vapor  pressure  of  air  (mb) 

C  ■  Cloud  amount  in  tenths 

h  **  Average  cloud  height  in  thousands  of  feet. 

To  expedite  the  use  of  Anderson's  equation,  Jcmes  [Ref.  3]  prepared  a 
nomograph  for  a  normal  range  of  values  (Fig.  2).  Cloud  height  and 
cover  are  expressed  in  terms  of  Lumb's  Cloud  Category.  Figure  (2)  was 
entered  using  sea  surface  temperature.  By  proceeding  horizontally  to 
the  atmospheric  vapor  pressure  and  moving  vertically  to  the  cloud  cover, 
then  horizontally  to  the  right,  the  value  of  effective  back  radiation 
was  found.  By  doing  this  for  several  temperatures  in  the  range  of 
observed  values  and  plotting  the  results  it  was  possible  to  linearize 
the  effective  back  radiation  as  a  function  of  sea  surface  temperature. 

D.  EVAPORATION 

Formulae  for  computing  evaporation  vary  from  simple  expressions 
relating  evaporation  to  vapor  pressure  differences  and. wind  speed  alone, 
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to  complex  considerations  involving  aerodynamic  surfaces,  the  occurrence 
of  spray,  and  the  vertical  profile  of  wind  and  vapor  pressure.  Laevastu 
[Ref.  14]  concluded  that  a  modification  of  the  formula  of  Rohwer  [Ref.  20] 
is  the  best  available  for  computation  of  evaporation  from  the  sea  surface. 
The  formula  as  modified  is: 

Qe  -  2.46  (0.26  +  0.04W)  (ew  -  e&)  (2-4) 

where  ew  *»  Saturated  vapor  pressure  at  sea  surface  temperature 
W  ■  Wind  speed  in  knots. 

The  case  where  eg  is  greater  than  ew  produces  condensation  and  a  heat 
gain  by  the  water.  However,  this  is  a  rather  rare  occurrence  in  mid 
latitudes.  The  above  formula  was  used  to  compute  Qe .  Since  ew  is  the 
variable,  it  was  necessary  to  use  another  of  James’  nomograms  (Fig.  3). 

By  entering  the  abscissa  with  sea  surface  temperature  (Tw)  and  rising 
vertically  to  the  diagonal  line,  e was  found  by  interpolating  between 
the  e  lines..  By  using  the  same  procedure  described  in  the  previous 
paragraph,  it  was  possible  to  linearize  ew  as  a  function  of  Tw  over  the 
limited  range  of  observed  values  of  sea  surface  temperature.  By  sub¬ 
stituting  this  linear  function  for  e into  the  above  formula,  the  energy 
loss  caused  by  evaporation  was  computed  as  a  function  of  sea  surface 
temperature . 

E.  SENSIBLE  HEAT 

Sensible  heat  transfer  Qh  is  a  function  of  the  vertical  temperature 
gradient  and  the  degree  of  turbulence  or  eddy  conductivity  above  the  sea 
surface.  Because  of  the  lack  of  information  on  the  eddy  conductivity,  it 
is  general  practice  to  relate  sensible  heat  transfer  to  evaporation. 
Assuming  that  evaporation  and  conduction  of  specific  heat  are  similar 
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processes,  Bowen  [Ref.  30]  derived  the  following  ratio  for  the  two 
processes  : 


Oh.  .  R  „  K  -(Ty  -  Ia)  . 

Qe  (e  -  e  )  1000, 

W  a 


(2-5) 


where  R  ”  Bowen  ratio 

Pa  ■  Atmospheric  pressure  (mb) 

K  ■  Constant  approximately  0.61  under  normal  conditions. 
Neglecting  the  pressure  term  which  has  only  a  small  effect,  and 
substituting  for  Qe  in  the  above  equation,  the  transfer  of  sensible 
heat  can  be  expressed  as : 

Qh  -  0.83  (0.26  +  0.04W)  (Tw  -  Ta).  (2-6) 

This  equation  is  applicable  where  colder  air  overlies  warmer  water. 

For  the  reverse  situation  stability  reduces  the  transfer  of  heat  and 
according  to  Laevastu  [Ref.  14]  the  appropriate  relation  is: 

Qh  -  0.036W  (Tw  -  Ta).  (2-7) 

The  cases  where  Ta  is  greater  than  Tw  are  not  numerous  and  in  light  of 
the  small  relative  magnitude  of  Qh ,  it  is  felt  that  equation  (2-6)  is 
a  good  approximation.  Since  this  equation  is  already  a  linear  function 
of  water  temperature,  no  further  manipulation  is  necessary.  James 
[Ref.  3]  has  plotted  a  nomograph  for  finding  latent  and  sensible  heat 
transfer  based  on  the  above  equations  (Fig.  4). 

The  sum  of  all  the  above  individual  linear  functions  of  sea  surface 
temperature  is  the  total  heat  transfer  across  the  air-sea  interface  and 
may  be  expressed  as  : 


Q  -  A  -  BT. 

The  subscript  has  been  dropped  since  only  water  temperatures  will  be 
discussed  from  this  point. 
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An  equilibrium  temperature,  Te ,  can  be  defined  such  that,  for  T  equal 
to  Te,  Q  is  equal  to  zero.  As  will  be  shown  later  in  the  development  of 
the  model,  a  ster.dy  state  condition  is  assumed.  For  steady  state 
conditions  to  exist,  the  ambient  sea  temperature  must  equal  the  equilib¬ 
rium  temperature.  Hence  the  constant  A  is  modified  such  that: 

A  -  BT“  "0  or  A  “  BT“ 
where  T“  •  Ambient  sea  temperature. 

This  in  effect  neglects  the  seasonal  heating  and  cooling  of  the  ocean. 
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III.  MODEL  DEVELOPMENT 


A.  GENERAL 

The  general  steady-state  heat  equation  for  an  advectivo  system  can 
be  written: 

“i  !r  -  se  C"i  Ir) ' < <T>  (3-1) 

i  1  /  i 

where  i"l,2,3  corresponds  to  the  x,y,z  direction  respectively  and 
repeated  indices  imply  summation,  represents  the  velocity  component. 
Dj,  is  the  eddy  diffusion  coefficient  for  heat.  Q(T)  accounts  for  any 
heat  gain  or  loss  by  the  system  and  may  be  expressed  as  a  function  of 
temperature.  This  allows  for  a  non-conservative  system  in  which  heat 
is  allot  to  be  lost  or  gained  from  the  body  of  the  fluid.  This 
equation  assumes  the  mean  heat  flux  is  proportional  to  the  negative 
temperature  gradient  (heat  flows  from  high  heat  to  low  heat)  where  the 
proportionality  factor  D^  is  the  diffusion  coefficient.  This  is  the 
so  called  Fickian  diffusion  where 

u. *T'  -  D.  §  (3-2) 

u^'  refers  to  the  turbulence  component  of  velocity,  T'  refers  to  the 
turbulent  temperature  fluctuation,  and  the  over  bar  signifies  a  time 
mean  average. 

B.  ASSUMPTIONS 

In  order  to  obtain  an  analytical  solution  for  the  heat  equation  it 
is  necessary  to  idealize  the  physical  situation  to  make  the  mathematics 
more  tractable.  The  following  simplifying  assumptions  will  be  made  in 
the  derivation  of  the  proposed  model. 
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1.  Only  a  uniform  mean  oceanic  longshore  current  will  be  con¬ 
sidered.  Additional  velocity  resulting  from  the  volume  discharge  of 
the  outfall  will  be  neglected.  This  is  reasonable  since  the  outfall 
only  affects  current  to  a  measurable  degree  in  the  immediate  vicinity 
of  the  boil.  The  velocity  is  also  assumed  uniform  over  depth  in  the 
diffusing  layer.  To  simplify  the  development,  the  x-axis  will  always 
conform  to  the  direction  of  flow,  although  the  results  could  easily  be 
extended  to  a  more  general  formulation  of  the  currents. 

2.  A  two  layered  system  separated  at  a  reference  depth  is 
considered.  There  is  complete  vertical  mixing  assumed  in  the.  top 
layer  and  no  mixing  between  the  two  layers.  This  closely  resembles 

the  thermal  structure  found  in  shallow  coastal  water  in  which  the  upper 
layer  is  rather  thoroughly  mixed  by  wind  and  wave  action  and  diurnal 
heating  and  cooling.  A  steep,  well-developed  thermocline  is  typical  of 
summer  conditions  in  the  temperate  zone  under  consideration.  The  water 
beneath  the  thermocline  tends  to  be  isothermal.  The  steep  thermocline 
tends  to  inhibit  vertical  mixing.  Wada  [Ref.  28]  found  that  the  hori¬ 
zontal  thermal  diffusivity  is  50  times  greater  than  the  vertical  thermal 
diffusivity .  Therefore  it  is  reasonable  to  assume  that  the  reference 
depth,  once  established,  can  be  assumed  to  be  a  constant  throughout  the 
horizontal  field. 

3.  As  already  mentioned,  the  ambient  sea  temperature  is  assumed 
equal  to  the  equilibrium  sea  temperature,  and  the  heat  budget  is  assumed 
a  linear  function  of  t -mperature . 

4.  A  steady-state  situation  exists. 

5.  The  horizontal  diffusion  coefficients  are  constants.  This 
assumption  greatly  simplifies  solving  of  the  differential  equation.  • 
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Although  Investigations  have  been  made  into  the  dependency  of  the  diffusion 
coefficients  on  distance  from  the  source  and  velocity  components,  none  of 
the  results  have  gained  universal  acceptance. 

6.  Ail  boundary  effects  except  for  the  surface  will  be  neglected. 
The  development  assumes  an  infinite  expanse  of  water  of  constant  depth. 

This  means  that  the  coastline,  shoaling  bottom  and  surf  zone  are  ignored. 

7.  The  outfall  will  be  treated  as  a  line  source  in  accordance 

with  the  two  dimensional  assumption.  This  is  a  fairly  realistic 
assumption  as  compared  to  the  actual  observations.  . 

Ui' 

C.  DEVELOPMENT 


Incorporating  the  above  assumptions,  the  governing  equation  reduces 
to  the  two-dimensional,  steady-state  equation  for  diffusion  of  heat  in 
a  moving  fluid: 


Dx 


a2T 

Sx2 


b2r 

-  Dy  — r 

dy2 


_ _ 9 _ 

D  C  Hw 
o  p 


(3-3) 


where  U(x)  »  Mean  current  (along  x-axis)  » 

Q  »  Net  gain  or  loss  of  heat  across  the  surface 
Dx  Diffusion  coefficient  in  x-direction 

I 

Dy  *  Diffusion  coefficient  in  y-direction  Ij  , 

oQ  *  Density  of  fluid,  a  constant 

c  *  Specific  heat  of  the  fluid 
P 

Hw  m  Specified  reference  depth. 

The  first  term  of  the  equation  represents  advective  heat  transport, 
the  next  two  terms,  the  thermal  diffusion  across  the  temperature 
gradients,  and  the  right  hand  term,  loss  of  heat  across  the  air-sea 
interface . 
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where  To  -  Outfall  (boil)  temperature 

T®  •  Ambient  temperature  ■  Equilibrium  temperature. 

Then 


and 


To  '  (T°  - 

■>*>£ 

S2T 

9 

-  (To 

-T-)& 

So 

So'2 

S2t 

o 

-  (To 

.  S2t 

-  T“)  — 7 

bri 

Sr] 

(3-8) 


Substituting  these  partial  derivatives  into  (3-6)  so  that 

>2T' 


U'  (To  -  T=°)  §-  *  <To 


Simplifying  this  equation, 


-  (if  ♦  if)  -  *•  -  »• 

W  Sri  ' 


T'  (To-T®)  +  T»l. 
(3-9) 


at'  /sv  , 

u  V' 


2t •  a'  -  B’T® 


-  B’T* 


(3-10) 


To  -  T® 

But  since  it  was  assumed  A  -  BT00  *  0  it  follows  that  A'  -  B'T®  **  0  and 

.  / 

the  equation  becomes 

*  ,%2_,  ,«*. 

(3-11) 


S1  -  <&■ +  iP ■  -8'1' 


Now  consider  the  transformation  from  a  stationary  to  a  moving  coordinate 
system.  Defining  another  new  variable 

?  -  a  -  U't  - - - - - - -  (3-12) 

where  t  is  an  arbitrary  time  variable. 

St 


Then  v1  -  I  and 
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0  so  that 
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This  transformed  equation  defines  a  new  temperature  field,  $  in  a  moving 

plane  in  which  the  advection  term  has  been  removed  and  4  is  being  diffused 

and  decayed  radially.  Hence,  it  is  convenient  to  work  in  cylindrical 

2  2  2  -1 

coordinates  with  4  (r,  6)  where  ?  +  r|  *  r  and  9  ■  tan  (n/§).  This 


transformation  is  immediately  written  down  as 

i±  +  l*i  +  i.  _  k2*  .  o 

3r2  ‘  S  r2  382 


(3-19) 


but  since  in  this  case  the  transformed  temperature  field  is  symmetrical 
about  r  ■  0,  the  function  4  must  satisfy  the  simpler  total  differential 


equation : 


d2?  .  1  d$  ,2 

— s-  +  —  - - k 

,2  r  dr 
dr 


0. 


(3-20) 


D.  SOLUTION 

The  differential  equation  above  is  recognized  as  a  Bessel  form 
having  the  general  solution  [Ref.  22] 

4  -  C1  JQ(kir)  +  C2  Yo(kir),  (3-21) 

or  in  terms  of  the  modified  Bessel  functions  and  re-expressed  in  T, 

T  -  (To  -  T»)  eD'?/2  [cL  IQ(kr)  +  C2  KQ(kr)]  +  T«.  (3-22) 

The  above  general  solution  must  satisfy  the  far  field  (asymptotic) 
boundary  condition, 

T - >  T=°  and  r  - >  ® .  (3-23) 

It  can  be  shown  that  of  the  two  functions  I  and  K  ,  only  K  will 

o  o^o 

satisfy  this  condition.  Accordingly,  ■  0  and 

T  -  (To  -  T»)  eU’?/2  [c  K  (kr)]  +  T®.  (3-24) 
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To  solve  for  C ^  the  rear  field  boundary  condition. 


-  DX  (2nr)  ^  ^ -  as  r  - >  0  (3-25) 

dr  Hwc  o 
P  ° 

is  applied  This  condition  expresses  the  fact  that  the  rate  at  which 
excess  heat  is  added  by  the  effluent,  qo,  is  equal  to  the  diffusive 
transport  of  heat  radially  outwards  near  the  axis  r  ■  0.  Here  it  is 

assumed  that  as  r  - >  0,  diffusion  is  constant  in  all  directions  arid 

Dx  represents  the  diffusion  coefficient.  Advection  is  neglected  very 
near  the  source  This  in  effect  neglects  the  conservation  of  mass  from 
the  volume  discharge,  and  leads  to  an  invalid  solution  in  the  immediate 
vicinity  of  the  outfall.  Although  this  appears  to  be  an  unrealistic 
boundary  condition,  in  practice  it  does  not  seem  to  be  any  great 
limitation.  By  applying  the  above  boundary  condition  it  is  found 
that ; 


c  .  _ 3°_ 

2  2ttHwDx  (To  -  T») 

and  the  solution  becomes: 


(3-26) 


T-e"'!/25FRk  Vkt>  +  T"-  «-27> 

’  .  y* 

Reversing  the  previous  simplifications  and  transformations  the 
solution  can  be  written: 


eu(x-ut)/2Dx  f _ aoj'N  ryuL. 

V2ttHwDx/  o  LvADx 


1/2 

B  >v  f(x-Ut) 

o  c  Hwy  \  Dx 
o  p 


2 


(3-28) 


The  above  solution  represents  a  steady-state,  non-conservative 
dispersion  model  for  heated  effluent  from  an  ocean  outfall  treated  as 
a  line  source.  The  behavior  of  the  solution  can  be  examined  by  con¬ 
sidering  the  components  of  the  solution.  The  exponential  term  increases 
in  the  downstream  direction  and  decreases  in  the  upstream  direction.  This 
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represents  the  influence  of  advection  and  is  the  term  in  which  the 
current  makes  its  largest  contribution  to  the  solution.  The  Bessel 
function  term  represents  the  radial  diffusion  of  heat  and  the  tempera¬ 
ture  field  decay.  The  combination  of  the  two  terms  plus  the  constant 
represent  the  difference  between  the  temperature  and  the  ambient 
temperature.  The  temperature  unrealistically  goes  to  infinity  at  the 
source,  but  the  solution  was  not  expected  to  be  valid  in  the  immediate 
vicinity  of  the  boil.  Por  downstream  values  of  x,  the  exponential  and 
Bessel  function  terms  vary  inverse ly  with  each  other;  the  Bessel  function 
term  decreases  at  a  slightly  greater  rate  thaft  the  exponential  term 
increases  resulting  in  a  decrease  in  temperature  with  distance.  In 
the  upstream  direction  as  the  distance  increases  both  terms  approach 
zero  and  the  temperature  approaches  the  ambient  temperature  more 
rapidly. 
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IV.  MODEL  APPLICATION 


The  model  was  calibrated  and  verified  using  the  data  from  the 
Huntington  Beach  power  plant.  The  temperature  and  current  surveys 
{Refs.  16  and  17]  were  conducted  at  the  different  tidal  conditions 
listed  below  during  the  late  summer  and  early  fall  of  1963  and 
16  February  1965  as  follows: 

26  July  1963  A  period  of  small  diurnal  and  small 

27  July  1963  semi-diurnal  tidal  components. 

16  August  1963  Periods  of  large  diurnal  and  large 

17  August  1963  semi-diurnal  tidal  components. 

16  February  1965  *. 

28  August  1963  A  period  of  large  diurnal  and  small 

semi-diurnal  tidal  components. 

4  October  1963  A  period  of  small  diurnal  and  large 

semi-diurnal  tidal  components. 


The  measurements  were  generally  made  only  during  daylight  hours. 
Temperature-measuring  instruments  used  on  these  surveys  were  a  con¬ 
ventional  bathythermograph  (MBT)  ,  an  electronic  bathythermograph  which 
is  capable  of  being  towed  at  fixed  depths  as  well  as  being  lowered 
vertically,  bucket  thermometers  and  an  airborne  radiation  thermometer. 

surements  were  made  by  placing  drogues  in  the  vicinity  of 
and  tracking  their  progress.  The  discharge  temperature, 
jlume,  and  meteorological  data  were  obtained  from  the 
lifornia  Edison  Company.  Additional  meteorological  data 


Current  me^ 
the  outfalJ 
discharge  v 
Southern  Ca 


were  obtained  from  the  City  of  Huntington  Beach  [Ref.  7]. 


An  evaluation  of  the  data  for  each  period  enabled  values  to  be 


found  for  To,  T=°,  Hw,  U,  qo,  and  B.  Values  for  T  were  presented  in 
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surface  temperature  contour  plots.  Ambient  sea  surface  temperature  (T=°) 
and  the  temperature  of  the  boil  (To)  were  taken  directly  from  the  surveys. 
The  value  for  the  current  (U)  was  based  on  an  analysis  of  the  drogue 
tracks.  These  tracks  showed  highly  variable  currents  in  most  cases. 
Complete  reversal  of  current  direction  within  a  period  of  two  to  three 
hours  was  not  uncommon.  The  currents  appeared  to  have  a  high  correlation 
with  the  tidal  cycle.  The  atmospheric  contribution  (B)  was  computed 
using  the  heat  budget,  and  the  excess  heat  added  by  the  outfall  (qo)  was 
estimated  by  taking  the  product  of  the  temperature  difference  (To  -  T=>) 
and  the  volume  discharge.  The  depth  of  the  temperature  field  (Hw)  was 
determined  from  the  actual  thermal  structure  as  measured  by  bathy¬ 
thermographs.  There  was  a  certain  amount  of  subjectiveness  in  the 
selection  of  a  single  Hw  for  each  period.  In  order  to  test  the  repeat¬ 
ability  of  selecting  Hw,  three  individuals  made  independent  observations. 
All  observations  were  within  20  percent.  To  solve  for  Dx  and  Dy,  it  was 
assumed  that  diffusion  is  horizontally  isotropic  (i.e,  Dx  **  Dy) .  This 
seems  reasonable  and  any  differentiation  in  the  diffusion  coefficient 
would  be  pure  speculation.  By  substituting  the  coordinates  and  tempera¬ 
ture  from  several  locations  on  each  surface  temperature  plot  into  the 
model  solution,  transcendental  equations  with  Dx  as  the  unknown  were 
established.  These  transcendental  equations  were  solved  and  the  value 
of  Dx  which  gave  a  "best  fit"  to  the  temperature  field  was  chosen  for 
each  survey.  By  using  an  exponential  approximation  for  large  values  of 
the  Bessel  function  argument,  a  closed  solution  for  Dx  was  possible. 

This  first  approximation  of  Dx  was  used  as  a  check  against  the  repre¬ 
sentative  value  chosen.  Substituting  Dx  back  into  the  model  solution, 
a  surface  temperature  field  was  generated  for  each  survey  period. 
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Comparisons  of  Che  generated  fields  and  the  actual  measured  fields 
may  be  seen  in  Figs.  6  through  12.  The  original  data  were  presented 


in  the  engineering  system 

of  units 

(i  .e  .  ,  feet ,  °F) .  For 

ease  in  com- 

parison,  the  g 

anerated  temperature 

fie  Ids  are 

presented  in 

these  same 

units.  Table 

I  lists  the 

parameters  used  in 

the  model  and 

!  the  values 

of  Dx  found . 

TABLE  I 

DATA 

SUMMARY 

- 

Dy 

4 

B  x  10 

qo 

Date 

Hw 

U 

O 

1  1 

2 

m  / sec 

m 

m/ min 

ca l/m“/hr 

cal/hr  x  10 

26  July 

.87 

3.00 

4.1 

3.92 

8.0 

27  July 

1.39 

3.60 

2.0 

3.06 

5.6 

16  August 

1.95 

3.00 

4.8 

3.06 

8.8 

17  August 

1.17 

2.20 

5.8 

.  3.06 

8.0 

28  August 

1.53 

3.60 

5.0 

2.86 

11.8 

4  October 

2.08 

3.00 

3.0 

2.86 

7.7 

16  February 

.78 

2.50 

5.0 

3.06 

4.0 

The  differences  between  the  observed  temperature  and  the  generated 

fields  can  be  attributed  mainly  to  irregular  current  patterns.  The 
model  is  constrained  to  a  steady  state  and  only  uniform  velocity  over 
the  entire  field  is  considered.  Internal  waves  and  upwelling  or  down- 
welling  cause  further  irregularities  it.  the  temperature  pattern.  Figure 
10  is  a  good  example  of  this  possibility.  In  all  of  the  figures,  the 
top  border  is  the  shoreward  direction.  The  effect  of  the  coast  and 
shoaling  may  be  seen  in  almost  all  of  the  figures.  A  further  cause  of 
differences  may  be  natural  oceanic  turbulence.  Roden  [Ref.  19]  con¬ 
cluded  that  random  temperature  fluctuations  of  ,5°F  may  be  expected. 
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Roden's  data,  actually  being  derived  from  deep  sea  investigations,  may 
in  fact  be  conservative  for  nearshore  regions  where  about  1°F  could  be 
a  significant  temperature  variation. 

A  further  comparison  was  made  to  determine  the  relative  importance 
of  the  heat  budget  contribution.  Using  the  data  for  27  July,  the  model 
was  applied  along  the  y-axis  of  the  field  both  including  and  excluding 
the  influence  of  the  heat  budget.  Figure  13  is  a  plot  of  the  temperature 
vs.  distance  along  the  y-axis  for  both  of  the  above  cases.  The  predicted 
temperatures  are  almost  identical  in  the  vicinity  of  the  boil  and  differ¬ 
ing  slightly  in  the  downstream  direction;  the  curve  representing  the 
temperatures  with  the  heat  budget  included  are  only  a  fraction  of  a 
degree  less  than  those  with  the  heat  budget  excluded.  This  is  because 
the  term  involving  B  in  the  Bessel  function  argument  is  several  orders 
of  magnitude  smaller  than  the  velocity  term  in  the  argument.  Thus,  at 
least  in  this  case,  the  heat  loss  term  could  be  neglected  without 
appreciable  error  in  the  results. 

The  effect  of  increasing  or  decreasing  the  velocity  and  the  dif¬ 
fusion  coefficient  was  examined.  It  was  found  that  the  areas  enclosed 
by  the  isotherms  varied  inversely  with  the  diffusion  coefficient.  The 
model  appeared  to  be  most  sensitive  to  changes  in  the  velocity  (U) .  A 
decrease  in  velocity  not  only  widened  the  field  as  would  be  expected, 
but  also  lengthened  it.  This  is  because  the  effect  of  the  current  is 
to  more  rapidly  disperse  the  heat,  decreasing  the  patch  areas.  From  the 
relative  importance  of  the  two  parameters  in  the  exponential  and  Bessel 
function  terms  of  the  solution,  one  would  expect  an  equal  and  opposite 
effect.  However  the  presence  of  Dx  in  the  constant  term  reverses  the 
effect  of  the  contribution  from  the  exponential  and  Bessel  function  terms. 
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The  effect  of  increasing  the  excess  heat  added  by  the  effluent  (qo) 
is  obviously  to  increase  the  areas  enclosed  by  the  isotherms.  Increasing 
the  reference  depth  (Hw)  would  have  the  opposite  effect.  The  contribution 
of  Hw  to  the  Bessel  function  term  is  negligible. 

An  examination  of  the  fields  generated  by  the  model  indicates  that 
near  the  boil  the  isotherms  are  more  circular  and  become  more  elongated 
as  the  distance  increases.  This  is  also  observed  in  the  surveyed 
temperature  field  and  lends  credence  to  the  second  boundary  condition 
applied  to  the  model  solution. 

An  attempt  was  made  to  correlate  the  parameters  used  in  the  appli¬ 
cation  of  the  mode  1,  part icularly  how  Dx  might  functionally  vary. 
Unfortunately  no  discernible  correlations  were  found.  It  is  felt  that 
this  is  not  a  result  of  discrepancies  in  the  model  but  reflects  the  lack 
of  a  sample  of  adequate  size.  The  lowest  diffusion  coefficient  was 
found  for  the  coldest  month  and  the  largest  diffusion  coefficient  was 
found  for  the  warmest  time  of  year  in  which  the  survey  was  conducted. 

This  might  indicate  that  Dx  is  a  function  of  the  thermal  structure  which 
is  related  to  Hw  and  qo.  However,  there  is  notable  exception  to  such  a 
correlation.  Another  low  diffusion  coefficient  was  found  for  the  26  July 
survey  during  which  time  the  winds  were  blowing  at  18  knots  and  white 
caps  were  being  generated.  This  was  the  most  severe  weather  condition 
of  all  the  surveys.  This  high  input  of  turbulent  energy  would  very 
definitely  modify  the  thermal  structure  which  in  turn  would  effect  Dx. 
Conversely,  the  sea  was  virtually  smooth  for  the  16  February  survey  for 
which  the  lowest  diffusion  coefficient  was  found.  More  data  will  have 
to  be  taken  and  analyzed  before  meaningful  correlations  can  be  established. 
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.  CONCLUSIONS  AND  RECOMMENDATIONS 


The  dispersion  model  developed  provides  a  method  of  predicting  the 

|  ‘ 

temperature  influence  from  a  heated  effluent  along  an  open  coast  using 
known  or  assumed  atmospheric  and  oceanic  parameters.  Diffusion  coeffici 
ents  were  calculated  in  the  calibration  of  the  model  as  compared  to  the 
field  data  at  iuntington  Beach.  The  value  of  Dx  varied  from  0.78  -  2.08 
m^/sec. 

An  attempt  was  made  to  correlate  the  diffusion  coefficient  to  the 
mean  velocity,  heat  input  and  thermal  structure.  No  correlation  was 
discernible;  this  is  attributed  to  the  lack  of  a  large  enough  sample. 

The  results  of  the  application  of  the  model  tend  to  over  predict  the 
temperatures  in  the  downstream  direction  and  under  predict  the  tempera¬ 
ture  in  the  upstream  direction.  This  is  evident  in  nearly  all  of  the 
figures  and  can  be  attributed  to  the  steady-state  constraint  on  the 
model.  Current  direction  fluctuations  tend  to  spread  heat  out  more 

evenly  in  both  directions  as  indicated  by  the  data.  The  contribution 

I 

from  the  heat  budget  in  this  application  was  shown  to  be  relatively 
insignificant  and  not  critical  to  the  results  of  the  model.  There 

I 

may  be  situations  where  this  term  may  be  significant,  therefore,  for 
generality,  it  is  recommended  that  the  term  be  retained. 

It  is  recommended  that  further  research  be  conducted  to  refine  the 
model  presented  in  this  study.  Some  possible  considerations  are: 

Addition  of  a  third  dimension,  coupling  jet  dilution  to  the  present 
model,  incorporating  the  current  induced  by  the  volume  discharge,  and 

introducing  time  variability.  Incorporating  any  of  the  above  parameters 

! 

j 

i 

.  i 

t 
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will  introduce  new  terms  to  the  differential  equation  and  make  an 
analytical  solution  unlikely.  More  detailed  and  accurate  surveys  will 
also  be  necessary  to  calibrate  these  more  complicated  models. 

A  further  recommended  study  would  be  the  correlation  of  boil 
temperature  to  condenser  discharge  temperature.  This  could  be  done 
by  taking  temperatures  of  the  boils  at  a  number  of  different  power 
stations  hopefully  with  a  wide  range  of  outfall  depths,  outfall  volumes 
and  condenser  discharge  temperatures.  An  analysis  of  the  study  could 
provide  an  empirical  formula  for  the  boil  temperature  as  a  function  of 
condenser  discharge  temperature,  outfall  depth,  outfall  volume  and 
outfall  diameter.  This  could  be  incorporated  in  the  above  model  and 
would  allow  an  effluent  field  prediction  starting  at  the  point  of 
discharge. 
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TABLE  n 

Data  for  Computing  Solar  Altitude 
(from  James) 


SUN'S  (°)  DURATION  OF  DAYLIGHT  (HRS) 

DATE  DECLINATION  20  30  40  50  60  Lat(N) 


January 

1 

-23 

11 

10 

9 

8 

6 

15 

-21 

11 

10 

10 

9 

7 

February 

1 

-17 

11 

11 

10 

9 

8 

15 

-13 

11 

11 

11 

10 

9 

March 

1 

-9 

12 

11 

11 

11 

10 

15 

-3 

12 

12 

12 

12 

12 

April 

1 

3 

12 

12 

13 

13 

13 

15 

9 

13 

13 

13 

14 

14 

May 

1 

14 

13 

13 

14 

15 

16 

15 

19 

13 

14 

14 

15 

17 

June 

1 

22 

13 

14 

15 

16 

18 

15 

23 

13 

14 

15 

16 

19 

July 

1 

23 

13 

14 

15 

16 

19 

15 

21 

13 

14 

15 

16 

18 

August 

1 

17 

13 

14 

14 

15 

17 

15 

13 

13 

13 

14 

14 

16 

September 

1 

8 

13 

13 

13 

14 

14 

15 

3 

12 

12 

13 

13 

13 

October 

1 

-3 

12 

12 

12 

12 

12 

15 

-9 

12 

11 

11 

11 

10 

November 

1 

-15 

11 

11 

11 

10 

9 

15 

-19 

11 

11 

io 

9 

8 

December 

1 

-22 

11 

10 

10 

9 

7 

. — . -  - 

15 

-23 

11 

10 

9 

8 

6 
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TABLE  m 

Lumb's  Cloud  Classifications  for  Computing  Insolation 

(From  James) 


CATEGORY 

DESCRIPTION 

1 

Virtually  clear  sky,  less  than  two  eights  coverage. 

2 

Well-broken  low  clouds  with  little  or  no  medium  or  high 
cloud  (three  to  five  eighths). 

3 

Six  to  eight  eighths  of  cirrus  (not  cirrostratus). 

4 

Thin  layers  of  altocumulus,  six  to  eight  eighths. 

5 

Veil  of  cirrostratus  over  whole  sky  with  up  to  four 
eighths  low  clouds. 

6 

Seven  or  eight  eighths  of  stratocumulus  without  rain, 
with  or  without  some  cumulus  and  little  or  no  medium 
clouds. 

7 

Thick  altostratus  six  to  eight  eighths,  with  or  without 
layers  of  stratocumulus  beneath  some  rain. 

8 

Thick  layers  of  stratus  and  stratocumulus,  overcast, 
including  drizzle. 

9 

Thick  layers  of  nimbostratus,  overcast,  also  include 
medium  clouds,  and  rain. 
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(from  James) 


VAPOR  PRESSURE  -  •„ 

8  10  14  18  24 


DEW  POINT  TEMPERATURE  -  Td 
SEA  SURFACE  TEMPERATURE  -  Tw 

FIGURE  3 

Vapor  Pressure  Nomograph 


(from  James) 
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WIND  -  KTS 


i 

i 

i 


Latent  Heat  of  Evaporation  and  Condensation  and  Sensible  Heat  Transfer 

(from  James) 


SCALE  :  FEET 


_ FIGURE;  6 _ 

SURFACE  TEMPERATURES  26  JULY  1963 

D««jB7ma/»ec 


48 


SURFACE  TEMPERATURES 


16  AUGUST  1963 


Di  -  1.95  m2t  tec 


5 


1 


SCALE  x  FEET 


«fc“ 


FIGURE:  10 


SURFACE  TEMPERATURES 
28  AUGUST  1963 
Ox  =1.53  m4/»*c 
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SURFACE  TEMPERATURE 
Dx-2.08  mZ/*ec 


4  OCTOBER  1963 
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SURFACE  TEMPERATURES  16  FEBRUARY  1965 

Dx=.78  ma/*«c 
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Significance  of  the  Heat  Budget 
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